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Abstract 

Differential cross sections for tlie reaction 7 p ^ n tt"^ have been measured witli the CEBAF 
Large Acceptance Spectrometer (CLAS) and a tagged photon beam with energies from 0.725 to 
2.875 GcV. Where available, the results obtained here compare well with previously published 
results for the reaction. Agreement with the SAID and MAID analyses is found below 1 GeV. The 
present set of cross sections has been incorporated into the SAID database, and exploratory fits 
have been made up to 2.7 GeV. Resonance couplings have been extracted and compared to previous 
determinations. With the addition of these cross sections to the world data set, significant changes 
have occurred in the high-energy behavior of the SAID cross-section predictions and amplitudes. 
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I. INTRODUCTION 



The photoproduction of mesons has played a crucial role in the search for resonances 
beyond those found through analyses of pion-nucleon elastic scattering data. Cross section 
structures seen in kaon and eta photoproduction [l| have been interpreted as candidates for 
so-called "missing" resonances, excitations that are predicted by QCD-inspired models 
but expected to couple weakly to the pion-nucleon channel. 

The photoproduction of pions, though less likely to detect states not seen in pion-nucleon 
studies, is the most well-developed of the meson-photoproduction programs, having an ex- 
tensive database for which many single- and multi-channel fits are available 3|. The photo- 
decay amplitudes for non-strange resonances have been determined almost exclusively from 
this reaction However, while cross section data exist, they are quite sparse above an inci- 
dent photon energy =1.7 GeV, and have generally come from untagged bremsstrahlung 
measurements. As a result, all photo-decay amplitudes for the higher A^* states have an 
inherent uncertainty beyond any model-dependence due to the background-resonance ex- 
traction process. While some theory-based model dependence is unavoidable, cross sections 
measured precisely using a tagged-photon beam, with incident photon energies covering the 
full resonance region, will provide tighter and more reliable constraints for future analyses 
of the properties of excited nucleons. 

In this paper, we report measurements of the unpolarized differential cross sections for 7r+ 
photoproduction on the proton for Ej from 0.725 to 2.875 GeV. As a first step to gauge their 
influence, we have included these new cross sections in a multipole fit to all available data 
covering the resonance region. This task is aided by the inclusion of tagged neutral pion cross 
sections recently measured 5|] that span a range in E^ from 0.675 to 2.875 GeV. We have 
obtained a revised set of multipole amplitudes and have extracted photo-decay couplings for 
those states that couple strongly to the pion-nucleon final state. Using the revised multipole 
analysis, we have generated predictions for further measurements of polarization observables 
that should soon become available. 

The paper is laid out in the following manner: We give a brief background of the ex- 
perimental parameters for this study in Sec. [Tll An overview of the method used is given 
in Sec. IIIII The uncertainty estimates for the cross sections obtained are given in Sec. llVi 
The experimental results are described in Sec. |Vl Various fits to the data are described in 
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Sec. IVl| and the underlying multipole amplitudes and resonance contributions are displayed 
and compared to previous determinations in Sec. IVIIi Finally, in Sec. IVlllt we provide a 
brief summary of the results of the study and consider what extensions of this work would 
be particularly helpful in the future. 



II. EXPERIMENT 



The differential cross sections for the reaction ^ p ^ n tt^ were measured with the 
CEBAF Large Acceptance Spectrometer (CLAS) 6| and the bremsstrahlung photon-tagging 
facility ( "photon tagger" ) [3] in Hall B of the Thomas Jefferson National Accelerator Facility 
(JLab) as part of a set of experiments running at the same time with the same experimental 
configuration (cryogenic target, tagger, and CLAS) called the "glc" run period. The cross 
sections were part of a program of meson photoproduction measurements undertaken using 



CLAS and the photon tagger 
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The data described here were obtained in sets of data runs with differing energies for 
the electron beam incident on the photon tagger. The two incident electron energies were 
2.445 and 3.115 GeV. Moreover, the 3.115 GeV data runs were taken with either the full 
photon-tagger focal plane ("3.115-fuU") or higher-photon-energy half of the photon-tagger 
focal plane ("3.115-half") in operation. Thus, for example, during the 3.115-half running, 
data was accumulated only for the higher-energy half of the available photon energies in 
order to increase statistics for data collected at those higher energies. 

The produced tagged photons impinged on an 18-cm-long liquid-hydrogen target placed 
at the center of CLAS. This target was enclosed by a scintillator array (called the "start 
counter," described in Ref. [l6]) that detected the passage of charged particles into CLAS 
from the target. The event trigger required the coincidence of a post-bremsstrahlung electron 
passing through the focal plane of the photon tagger and at least one charged particle 
detected in CLAS and the start counter. Tracking of the charged particles through the 
magnetic field within CLAS by drift chambers |l3] provided determination of their charge, 
momentum, and scattering angle. This information, together with the particle velocity 
measured by the time-of-fiight system ISj and start counter, provided particle identification 
for each particle detected in CLAS and its corresponding momentum four-vector. 

The methods used for extracting the differential cross sections for 7?+ photoproduction 
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are presented in the next several sections. The technique is outhned initially, and then each 
step is described in further detail, with a summary provided of the data and tests that 
support the validity of the approach taken. 

III. DATA REDUCTION 

The technique for this analysis is very similar to that used previously in the analysis of 
the CLAS glc running period data for the reaction jp pir^ [5]. In that analysis, the 
recoiling proton from the photoproduction process was detected in CLAS and, assuming the 
two-body reaction 7 p — > pX (where X was the undetected tt"), yields were determined in 
the missing mass spectra for the reconstructed 7r°. 

In this analysis, similarly, the photoproduced 7r+ was detected in CLAS. Again assuming 
the two-body reaction 7 p ^ tt+X, where in the present case X was the undetected neutron, 
yields were determined in the missing mass spectra for the reconstructed neutron. However, 
while both the proton and tt~^ are positively charged particles, the CLAS detector response 
to the recoiling pions and protons was different (for example, the amount of energy deposited 
in the scintillators within the detector), which necessitated appropriate modifications to the 
previous analysis. 

For the data described in this paper, yields for the neutron were determined using the 
following steps: 

1. Identify the 7r+ in CLAS, determining the scattering angle and momentum. 

2. Sort the events in the resulting missing mass spectra into kinematic bins in 
incident photon energy and scattering angle 9'^^ , where 9'^^ is the center- 
of-mass angle of the 7i~^. 

3. Identify the missing mass peak for the neutron in each kinematic bin. 

4. Determine the yield for the neutron in each kinematic bin by subtracting the 
background beneath the peak. 

5. Correct the meson yield in each kinematic bin for spectrometer acceptance using 
a Monte Carlo simulation of the spectrometer acceptance. 
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FIG. 1: Particle identification spectrum obtained with CLAS, showing identifications provided by 
the GPID algorithm (discussed in the text) for all charged particles. 

6. Normalize the measured yield in each kinematic bin using a measured absolute 
photon flux normalization procedure, thereby determining the differential cross 
section for that bin. 

In the following sections, each of these steps is described. Also presented are sample 
results, and, in some cases, tests that establish the validity of the procedures used. 



A. Particle identification and kinematic variables 



The tracking information provided by the drift chambers within CLAS gave momentum 
and scattering angle information on charged particles scattered within the detector volume. 
Time-of-flight and start counter information, coupled with the track information provided 
by the drift chambers, determined particle velocity and momentum. 

Particle identification in this analysis was performed using the GPID algorithm 19(]. The 
method uses the momentum of the detected particle, and sequentially calculates trial values 
of the velocity (3 for all possible particle identities. Each one of the possible identities is 
tested by comparing the trial value of (3 for a given particle type to the empirically measured 
value of (3 (as determined by CLAS tracking and time-of-flight information). The particle is 
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assigned the identity that provides the closest trial value of (3 to the empirically measured 
value of [3. For example, if the curvature indicates a positive particle, the [3 is calculated 
for j9, TT"*", and . Figure [T] shows the mass distribution of the identified charged particles. 
The GPID algorithm also attempts to find a matching photon in the tagging system for 
every charged particle detected in CLAS. A matched photon means that there was one and 
only one tagged photon in the trigger window, which, in this analysis, was 18 ns. Particles 
that were determined not to have a matching photon are considered to be a measure of the 
accidentals (to be described in more detail in the next subsection). 

CLAS is divided into six sectors in azimuthal angle. Geometrical fiducial cuts in each of 
the six sectors of CLAS were imposed on all pions. The region selected for accepting pions 
in each sector corresponded to a region of relatively uniform detection efficiency (constant 
to ± 3%) versus azimuthal angle. 

B. Missing mass reconstruction 

The momentum for the 7r+ was determined by the drift chamber system. The momentum 
determined by CLAS was corrected for energy loss in both the target cell and the start 
counter [2(3]. The scattering angle and momentum was used to calculate the missing mass 
based on the assumption that the reaction observed is vr+X. Based on this assumption, the 
missing mass spectrum in the full spectrometer acceptance for all photon energies is shown 
in Fig. [21 The neutron peak is clearly seen. 

Taking each vr"*" event that did not have a matching incident photon as noted above, and 
integrating over all of the out-of-time (not within the trigger coincidence window) incident 
photons for that event, determined the distribution of accidental coincidences between CLAS 
and the photon tagger. This assumes that coupling the out-of-time tagger hits to unmatched 
pion created a fair representation of the accidental coincidences between CLAS and tagger. 

C. Distribution of events into kinematic bins 

The events from both the 2.445 and 3.115 GeV data sets, constituting the full missing 
mass spectrum described in the previous section, were sorted into bins in incident photon 
energy and cos . The widths of these "kinematic bins" ( Ai?^ = 50 MeV in photon 
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FIG. 2: Missing mass spectrum obtained from the glc data set using CLAS, assuming the reaction 
7 p — > tt'^X . 

energy and Acos^Jj„ = 0.1) were chosen such that, in general, there were at least 1000 7r+n 
events in each kinematic bin. 

D. Neutron yield 

For each kinematic bin, the neutron yield was extracted by removing the background 
under the peak. We have proceeded with the assumption that the background in the missing 
mass spectra arises from two particular types of events: 

1. Events arising from accidental coincidences between CLAS and the photon tag- 
ger, as discussed in the preceding subsection. 

2. Events arising from two-pion photoproduction via the reaction 7 p — > n'^X, 
where X = pn~ 01 X = n7r°. 

The spectrum for accidental coincidences is determined by looking at events that fell 
outside the designated trigger window. To determine the two-pion background, data for the 
reaction 7 p — > p vr" 7r+ were selected by requiring that each particle in the final state had 
to be identified through normal particle identification procedures, that the same incident 
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FIG. 3: Neutron yield extraction for = 1.475 GeV. The background is represented as the 
shaded region. 

photon was chosen for each particle, and that the missing mass was consistent with zero. 
These selected data were used to determine the shape of the X = p7i~ and X = rni^ 
components of the background from two-pion photoproduction (due to A(1232) dominance, 
the contribution from the X = nir^ reaction was assumed to have the same shape as the 
X = p-K^ contribution). This shape was used to generate the background beneath the 
neutron peak, which was then subtracted from the neutron yield for each kinematic bin. 
The fractional uncertainty in the background beneath the peak was statistically added in 
quadrature to the uncertainty in the yield for each kinematic bin. In most cases {^93 %), the 
peak-to-background ratio was greater than 5 to 1; in all cases, the signal-to-background ratio 
was greater than 1.4 to 1. Fig. [3] shows an example of this background removal procedure 
for all kinematic bins with photon energy = 1.475 GeV. 



E. Acceptance and efficiency 

The spectrometer acceptance for charged pions was determined from the results of Monte 
Carlo simulations of the CLAS detector response to positive pions. As a preliminary test 
of the quality of the Monte Carlo representation of the CLAS response to 7r+, simulated 
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acceptances for tt"*" were compared to empirical measurements based on the reaction 7p 
p7r~7r"^ for most of the kinematic bins in this study (the empirical acceptance method is not 
useful for some regions of phase space due to limited statistics for those kinematic bins). 
Such an empirical check is practical for much of the phase space covered in this experiment 
due to the large number of events for that final state, and that all of the final products leave 
charged tracks in CLAS, making them easily observable. For the empirical comparison, in 
addition to the tt"^, the proton and 7r~ were required to be detected in the event and both 
were assigned the same photon. The same fiducial cuts applied to the tt"^ noted above were 
applied to both reconstructed and CLAS- identified tt^. A missing mass reconstruction from 
the kinematic information of the proton and vr^ was performed to determine if a 7r+ should 
have been seen in CLAS. The background beneath this peak was removed by subtracting a 
polynomial fit (order 3) to the background region from the spectrum. 

A comparison of Monte-Carlo-simulated events to actual data for the 7 p — > p 7r~ tt"^ 
reaction (re-binned as if the tt"^ came from the 7 p — > n tt"^ reaction channel) was performed. 
Simulated events were obtained by generating 10^ 7 p — >^ p 7r~ 7r+ events that were isotropic 
in phase space and then processed through a GEANT simulation of CLAS created by the 
Jefferson Laboratory GSIM working group. In addition to simulating the detector response, 
the GEANT simulation also included the effects of pion decay. In those kinematic bins 
where the acceptance was less than 10%, agreement between the empirical and Monte- 
Carlo-simulated acceptances was poor. Thus, an acceptance cut was applied such that only 
kinematic bins that had acceptances greater than 10%, and had no neighboring bins with 
acceptances less than 10%, were kept. In addition to this "10% criterion," the bins at 
cos^^jj^ > 0.9 and cosO^j^ < —0.9 were removed, since some portion of these bins would 
have had acceptances of zero due to the geometry of CLAS. The fraction of all kinematic 
bins rejected by the "10% criterion" was 0.195. 

The empirically-measured and Monte-Carlo-simulated acceptances agreed well when these 
conditions were applied. To quantify this agreement, an "acceptance ratio" was determined, 
defined as the ratio of the empirical acceptance to the Monte Carlo simulated acceptance for 
each photon energy and cos^^^^^ bin (with the acceptance cut applied). These acceptance 
ratios were placed in a histogram, and then fit with a Gaussian. The center of the Gaussian 
was 0.9997 and the standard deviation was 0.040, which affirms the validity of the Monte 
Carlo simulation of the response of CLAS to 7r+. 
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In addition to examining the ratio of the empirical acceptance to the Monte Carlo accep- 
tance, a standardized Gaussian distribution Zij was created by forming, for each kinematic 
bin, the difference of the Monte Carlo simulated acceptance eMc and the empirical accep- 
tance eE, with that difference divided by the combined acceptance uncertainty thusly: 

\"^3)mC+E 

where 



histogramed for each energy i and cos^^j^ j kinematic bin. These points are assumed to 
obey Gaussian statistics with a variance of one and a centroid located at exactly zero. If 
the centroid of the distribution has been "pulled" away from zero, that suggests the Monte 
Carlo acceptance results {Uj)Mc ^'^ approximate the empirical acceptance exactly. If 
the variance of the Zij distribution is less than one, then the uncertainties {(^ij)Mc+E ^'^^ 
large. Conversely, if the variance of the Zij distribution is greater than one, this suggests 
the uncertainties {.(^ij)Mc+E underestimated. 

The uncertainties of the Monte Carlo acceptance are assumed to be well-represented by 
the uncertainty appropriate for a binomial distribution: 



i^ijjMc = \ — ' 

V y-^^'tj) Thrown 

where {Nij)rpf^^^^^ is the number of events thrown in the ij kinematic bin. 

The mean of the standardized Gaussian distribution Zij was from equation ([1]) nearly 
equal to zero within uncertainties, 0.10 ± 0.09. The value of Xreduced ^^e Gaussian fit 
to the distribution, Xreduced ~ 0.86, was also reasonable. However, the standard deviation, 
1.29 ± 0.09, was larger than the optimal value of one, suggesting that the uncertainties 
o'mc+e were too small. When an additional 2% uncertainty was added in quadrature to the 
Monte Carlo uncertainty, the centroid, standard deviation, and Xreduced were 0.09 ± 0.07, 
1.02 ± 0.05, and 0.514, respectively. 

To test how far the Monte Carlo results were from optimal, we added 0.1% to the Monte 
Carlo efficiency. With this 0.1% shift to the Monte Carlo the centroid, standard deviation, 
and Xreduced wcrc (0.05±0.06, 1.02±0.05, 0.465, respectively). Since the values are consistent 
with the optimal values, we assume henceforth that the Monte Carlo acceptances agree very 
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well with the empirical acceptances when 2% additional uncertainty is added to the Monte 
Carlo acceptances. The remainder of this analysis assumes that it is appropriate to add 
this extra 2.0% uncertainty in quadrature to the Monte Carlo uncertainties on a bin-by-bin 
basis, and that has been done for each kinematic bin. 

Having confirmed the validity of the Monte Carlo representation of the CLAS response 
to vr"*", the acceptance results for the reaction 7 p ^ n Tr"*" were obtained by generating 
10^ events (weighted by the cross sections given by the SAID solution [5|]) that were then 
processed in the same manner as the 7 p ^ p vr^ 7r+ comparison reaction. These simulated 
acceptances were used to determine the differential cross sections reported here. 

F. Sector-by-sector comparison 

A sector-by-sector comparison of the differential cross sections was performed to check 
the consistency of the extracted cross sections. CLAS is constructed from six sectors which, 
ideally, should be identical. However, operationally, the response of each sector is different 
owing to various hardware circumstances, problems, and differences. The simulations de- 
scribed in the previous section incorporate knowledge of the various differences in the sectors 
in order to properly reproduce the CLAS response for each particle type. Since the Monte 
Carlo simulations should reflect sector-by-sector changes in the detector arising from, for 
example, holes in the drift chamber system due to broken wires and bad time-of-flight pad- 
dles, a sector-by-sector comparison of the differential cross sections inferred from the data 
obtained explores the reliability of the Monte Carlo with respect to these detector irregulari- 
ties. The results of this comparison indicated that variations attributable to sector-by-sector 
variations were less than 0.4 %, and much smaller than the uncertainty in the cross sections, 
thus confirming the validity of the simulated sector-by-sector response. 

A standardized Gaussian distribution for the sector-by-sector comparison was created by 
forming, for each photon energy, cos 6'^^ , and sector bin, the difference of the differential 
cross section in each sector to the sector average and dividing the result by the uncertainty. 

The resulting centroid, standard deviation, and xleduced standardized Gaussian 

distribution were 0.047 ± 0.021, 0.979 ± 0.018, and 1.01, respectively. Thus, while the 
xleduced Standard deviation of the Gaussian are reasonable, the centroid is somewhat 
smaller than the optimal value of zero. 
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To roughly estimate how far off the cross sections might be from the desired value for the 
centroid, we shifted the sector average by a factor of 0.996. The resulting modified centroid, 
standard deviation, and Xreduced^ were found to be, 0.003 ± 0.021, 0.985 ± 0.018, and 0.979, 
respectively. Since this small shift of 0.4 % in the sector average, a shift much smaller than 
the uncertainty for the cross section, produces parameters for the standardized Gaussian 
that are within optimal values, the non-shifted parameters are acceptably close to optimal. 

G. Bin migration 

To estimate the systematic uncertainty associated with bin migration, the acceptance 
and efficiency results calculated using SAID-weighted events were compared to acceptance 
and efficiency results using non-weighted events. Since the amount of the correction was 
found to be typically less than 2%, the systematic uncertainty associated with bin migration 
was assumed to be ignorable. 

H. Trigger inefficiency 

The determination of a charged particle trigger inefficiency for the glc data was performed 
by looking at data from a running period just preceding the glc period, the g2a running 
period. (The g2a running period is more fully described in Ref. [21,]). This running period 
had, in addition to the charged particle trigger, a photon trigger. The photon trigger required 
that there was a hit in any two sectors of the electro-calorimeters located downstream of 
CLAS in coincidence with a hit in the photon-tagger. By looking at g2a events that had 
a photon trigger and no charged trigger, yet had a vr"*" in the event, the inefficiency of the 
charged particle trigger in CLAS for 7r+ was determined. This correction was applied to 
each kinematic bin, and was always less than 1.0%. 

I. Normalization 

The absolute photon flux for the entire tagger photon energy range was determined by 
measuring the rate of scattered electrons detected in each counter of the focal plane of the 
bremsstrahlung photon tagger by sampling focal plane hits not in coincidence with CLAS. 
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The detection rate for the post-bremsstrahlung scattered electrons was integrated over the 
hfetime of the experiment and converted to the corresponding total number of photons on 
target for each counter of the tagger focal plane. The tagging efficiency was measured in 
dedicated runs with a total absorption counter (TAG) downstream of the cryogenic target, 
which directly counted all photons in the beam. The details of the method can be found in 
Ref. j22]. 

IV. UNCERTAINTIES 

We summarize here the various uncertainties present in the cross sections obtained in 
this work. 

• An overall estimated systematic uncertainty of 1% is taken as a very conservative 
estimate of all sources of trigger inefficiency, as described in section IllIHI 

• The uncertainties associated with the detector response, bin migration, and track re- 
construction are contained within the uncertainties associated with the Monte Carlo 
acceptance estimates described in subsection llll El These uncertainties are taken into 
account on a bin-by-bin basis. 

• The uncertainties associated with the background subtraction described in subsec- 
tion IIII Dl are purely statistical, and these were taken into account on a bin-by-bin 
basis. 

• The largest source of uncertainty in the photon flux normalization arises from the 
uncertainty in the measurement of the "tagger efficiency" jj], essentially a measure 
of the photon beam coUimation taken during normalization runs. The value of this 
tagger efficiency is dependent upon the positioning of the electron beam supphed by 
the accelerator on the radiator of the photon tagger, and will vary on a run-by-run 
basis determined by the run-by-run condition of the electron beam tune. With the 



procedure used to obtain the photon flux normalization [22| , the statistical uncertain- 
ties associated with the photon flux normalization are always far below 1% and, when 
considered with other uncertainties in the absolute normalization, are negligable. 
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• The systematic uncertainty of the absolute normalization is comprised of six parts; 
three of them do not vary over the running period, while the remaining three do. The 
following quantities vary over the running period: 

1. run-to-run variations in the normalized neutron yield unaccounted for by 
statistical uncertainties alone; 



2. uncertainty in the target density [23|; 

3. statistical uncertainty of the photon flux normalization. 

Table |T] shows contributions to the systematic uncertainties for quantities that varied 
over the running period. 

The following systematic uncertainties do not change over the running period: 



1. uncertainty in the liquid-hydrogen target-cell length, which was 0.3% 2J]; 

2. uncertainty associated with the tagger energy calibration (described in sub- 
section [Till]), which was less than 1%; 

3. uncertainty in the trigger inefficiency correction, which was less than 1%; 

After adding all of the systematic uncertainties in quadrature, the systematic uncertainty 
for the absolute normalization is 1%, 2%, and 4% for the 2.445 GeV, 3.115 GeV (full) and 
3.115 GeV (half) data sets, respectively. However, since combinations of more than one of 
these data sets was used to obtain the differential cross section for each kinematic bin, a 4% 
absolute normalization uncertainty is assumed for simplicity. 

TABLE I: Systematic uncertainties in the absolute normalization for quantities that varied over 
the running period. (The data set descriptions are discussed in section II.) 



Data set 


Run-to-Run 


Target density 


2.445 


0.9% 


0.1% 


3.115-full 


1.9% 


0.3% 


3.115-half 


3.6% 


0.3% 
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V. RESULTS 
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stained in this experiment are compared to the world 
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44 



45| 



m 



Figs, m — El along with a number of representative fits described below. The differential 



cross sections reported here are the first tagged TT~^n measurements above 780 MeV 45|. 



The cross sections are available in electronic form in Ref. 47|- The database entries in- 
clude the differential cross sections, as well as uncertainties (excluding the overall absolute 
normalization uncertainty), for each incident photon energy and cosO'^^^ bin shown in this 
paper. 

For a specific example of agreement with previous measurements, in Fig. H] we compare 
differential cross sections obtained here with those from the A2 collaboration of the MAMI- 
B group [45], at an energy common to both experiments. The CLAS data and the results 
from MAMI-B appear to agree well at this energy. 
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FIG. 4: The differential cross section for -fp — > Tr~^n at = 725 MeV versus pion center-of-mass 
scattering angle. Solid (dotted) lines correspond to the SAID FA08 (FA07) solution. Dashed lines 



give the MAID07 [46[ predictions. Experimental 



MAMI-B measurements (triangles ) at 723 MeV 



circles) are from Refs. 



30, 



3^, 



33, 



3S 



data are from the current (filled circles) and recent 



451 ]. Previous bremsstrahlung measurements (open 



43l ]. The data have been selected from energy bins spanning 



at most 3 MeV. Plotted uncertainties are statistical. 

More generally, as can be seen in Figs. S] — El agreement with previous measurements 
is good overall. The largest deviations generally occur at forward angles. Thus further 
measurements at more forward angles would be useful. While agreement with previous 
measurements is generealy good, even so, the data here extend measurements to higher en- 
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ergies with more complete angular coverage than obtained in those previous measurements. 



VI. AMPLITUDE ANALYSIS OF DATA 

We have included the new cross sections from this experiment in a number of multipole 
analyses covering incident photon energies up to 2.7 GeV using the full SAID database in 
order to gauge the influence of the present measurements, as well as their compatibility 
with previous measurements. A "forced" fit, which included the present dataset weighted 
by an arbitrary factor of 4, was compared to a standard fit. (The standard fit with normal 
weighting is called henceforth FA08.) The results with two different weightings were in good 
agreement, despite the CLAS data having a larger weighting. This agreement is not surpris- 
ing concidering the agreement of these new data with previously published measurements 
and that an older fit (FA07) was able to give a reasonable prediction for the previously 



published cross sections at all but the highest energies [48 1. 



In Table HTl we compare FA08 with two previous SAID fits (FA07 and FA06 [Sj) and also 



0: 



with the Mainz fit MAID07 j46| up to its stated center-of-mass energy W limit of 2 GeV 
{E^ = 1.65 GeV). The FA07 fit included LEPS Collaboration ir^p measurements [49]. These 
three solutions are compared with the data in Figs. [5] and [61 While the FA07 and FA08 
SAID fits agree well over the energy range of the Mainz fit, disagreements between the SAID 
and MAID fits are most pronounced at angles more forward than the CLAS data. Near its 
upper energy limit, the MAID07 solution also exhibits structure not seen in the data. 

Above 2.4 GeV, the new CLAS data reported here begin to depart from the FA07 pre- 
dictions. As a result, the new data presented here have resulted in adjustments of a number 
of parameters in the FA08 solution so that the new solution better reproduces the measured 
cross sections, which are significantly lower than the predictions given by FA07. 

In fitting the data, the stated experimental systematic uncertainties have been used as 



an overall normalization adjustment factor for the angular distributions [SOj. FA07 included 
all previously published data used in FA06 jsl, plus recent 7r°p differential cross sections and 



beam asymmetry S data from the LEPS Collaboration [49|. The MAID07 analysis does 
not include the recent tt^p measurements from CLAS [5] and LEPS [49], and has a center- 
of-mass energy limit of W = 2 GeV {E^ = 1.65 GeV). Presently, the pion photoproduction 
database below E^ = 2.7 GeV consists of 25639 data points that have been fit in the FA08 
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solution with = 54161. The contribution to the total iii the FA08 analysis of the 561 
new CLAS vr+n data points (e.g., those data points up to = 2.7 GeV) is 1407. 

TABLE II: comparison of fits to pion photoproduction data up to 2.7 GeV. Results are shown 
for three different SAID solutions (FA08, FA07, and FA06) recent MAID07. See text for details. 
Comparison includes all previous plus new CLAS vr+n measurements. 



Solution 


Energy limit 
(MeV) 


X^/Data 


Data 


FA08 


2700 


2.11 


25639 


FA07 


2700 


2.02 


24376 


FA06 


3000 


2.15 


25252 


MAID07 


1650 


7.38 


22621 



Multipoles from the FA08 fit are compared to the earlier MAID07 determinations in 
Figs. [7] and [HI Both FA07 and FA08 are quite similar, but significant differences between 
SAID and MAID in magnitude (e.g., i?2^^,M|^^, and E^!^) and W dependance (e.g., mI/^^, 
and MfL^) are seen. Given that large differences are not seen in the differential cross sections, 
further measurements of spin observables will be needed to better constrain these amplitudes. 
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FIG. 5: The differential cross section for — > n^n below = 2.7 GeV versus pion center-of- 



mass scattering angle. So 
lines give the MAID07 



id (dotted) lines correspond to the SAID FA08 (FA07) solution. Dashed 



461 ] predictions. Experimental data are from the current (filled circles) and 



previous measurements (open circles). The plotted points from previously published experimental 
data are those data points within 3 MeV of the photon energy indicated on each panel. Plotted 
uncertainties are statistical. 




FIG. 6: Fixed angle excitation functions for 7p — > n'^n. The pion center-of-mass scattering angle 
is shown. Notation as in Fig. [5j The plotted points from previously published experimental data 
are those data points within 2 degrees of the angle indicated on each panel. 
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FIG. 7: Multipole amplitudes from threshold to E.^ = 2.7 GeV for isospin 1/2. Solid (dashed) lines 
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show full r and partial widths for TttAt associated with the SAID vrA^ solution SP06 j5ll |. 
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With the addition of CLAS vr^p and Tr+n cross sections, the SAID solution at higher 
energies is now far more rehable than in previously published analyses. Based on the earlier 
SAID SM05 solution, the authors of Ref. 52| previously noted how well the single-pion 



component of the Gerasimov-Drell-Hearn (GDH) sum rule integrand reproduced the full 
result (including multi-pion and other-meson production). In Fig. [HI we extend this same 
comparison significantly beyond the 2 GeV range of the SM05 solution. As seen in the 
figure, the FA08 solution now agrees well with the MAID07 result, but extends that result 
to much higher E^. General agreement with the existing GDH data 53| is good. 
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FIG. 9: Single-pion photoproduction contributions to the 
from the SAID current (solid), recently published SM05 



Droton GDH sum rule Ao"3i 



analyses. GDH data from Ref. 
quadratures. 



52] (dotted), and MAID07 



0"3/2-0"l/2 



461] (dashed) 



531 ] . Plotted uncertainties are statistical and systematical added in 



For completeness, we provide in Fig. (TD] a comparison between the predictions for the 



beam asymmetry E from the FA07, MAID07, and 



AOS analyses and the experimental data 



for that variable from GRAAL [54], from DNPL [55], and from CEA 56] for the n^n 
reaction under study here. The agreement with the GRAAL data for S at 1.3 GeV is very 
good for both SAID solutions, while there are discrepancies at center-of-mass scattering 
angles greater than 75° between those data and the MAID07 predictions. All three analyses 
are seen to match the single S data point from CEA at 1.6 GeV, and both the FA07 and 
FA08 analyses provide reasonably good predictions for the DNPL data for S for positive 
pions at 2.1 GeV 55|, although the agreement is poorer for center-of-mass scattering angles 
greater than 75°. However, the data for S remain relatively sparse compared to the existing 
data for the differential cross sections. New data for S will help firm up the experimental 
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situation for this energy region, and a number of experiments are underway at Jefferson Lab 



to obtain such data for pions and other mesons 
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FIG. 10: Beam asymmetry S for 7p n^n at E-^ = 1300, 1600, and 2100 MeV vs. center-of-mass 
scattering angle. Solid (dotted) lines correspond to the SAID FA08 (FA07) solution. Dashed lines 

redictions. Experimental data are from GRAAL (filled circles) [s^, from 
5^ , and from CEA (open triangles) [56:]. Plotted uncertainties are statistical. 



give the MAID07 [46 
DNPL (open circles) 



Systematic uncertainties are taken into account in the fit (see text). 
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VII. RESONANCE COUPLINGS 

As in Ref. [sl, we have extracted resonance couplings from the modified fit (FA08) using 
a simple resonance plus background assumption, a form similar to that used in the MAID 
analysis, 

B{W) {l + tT^N)+TBw (4) 

where T^^jsi is the associated full pion-nucleon T-matrix and Tbw is a Breit-Wigner 
parametrization of the resonance contribution. With this model, resonance contributions 
have been determined and are listed in Table UTTl Values for the resonance mass Wr, width 
r, and branching fraction F^^f/r for the various resonances were taken from a recent SAID 
analysis of pion-nucleon elastic scattering data 5l|] • These couplings were also calculated in 
Ref. [5[ after the addition of vr^p photoproduction data reported in that reference. 

The function B{W) was fit to data over an energy range spanning the resonance position. 
In the MAID determination, B{W) was given by the Born term. Differences between the 
couplings quoted here and in MAID therefore reflect both the impact of the present data 
set and a model-dependent uncertainty associated with the resonance extraction procedure. 

Results based on a fit not including the present data set, presented in Ref. j^, generally fall 
within one to three standard deviations of the present values. This stability is to be expected; 
larger deviations may occur with the addition of forthcoming polarization measurements. 

However, the range of couplings given in Table IIIII requires further comment. The two 
resonances coupled to a -kN Su state are given very different estimates in the present analysis 
than those provided by the 2007 MAID fit and the PDG. The PDG range for the iV(1535) 
accounts for the large discrepancy that once existed between determinations based on vrA^ 
and rjN photoproduction fits. Whereas the present vrA^ estimate, the PDG central value, 
and older rjN photoproduction analyses agree on a value close to 100 GeV"^'^ x 10~^, the 
MAID 2007 value has now dropped to a value consistent with the 1996 SAID value js^]- 
This low value was criticized in a number of papers analyzing rjN photoproduction data 
measured at MAMI-B in Mainz 



From the plots in Figs. [7] and [HI a significant difference between the SAID and MAID 
fits exists in multipoles coupled to the vrA^ and Di^ resonances. This, combined with 
differences in the assumed background contribution, likely accounts for the variations seen 
in Table HITl Differences in the A^(1650) couplings are largely due to difficulties in separating 
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two nearby resonances in a single multipole. The present A^(1650) photo-decay amphtude 
is consistent with that found in Ref. j^, given the large errors. The statistical significance 
of any inconsistencies with the MAID analysis cannot be determined, as they have not 
presented any uncertainties for their estimates. 

Both the SAID and MAID values for the A^(1720) coupling are very different from the 
Particle Data Group (PDG) average. The PDG range does not even determine a sign for 
this coupling. As this state has the lowest vrA^ branching fraction listed in Table IIIIl a 
better determination may require a more favorable reaction or additional information on 
spin observables. Finally, we note that while the present SAID fit, the fit in Ref. 5|, and 
the PDG estimate for the A(1700) photo-decay amplitudes have remained relatively stable, 
the MAID 2007 value for Ai/a amplitude has nearly doubled the MAID 2003 result. This 
change has resulted in both the helicity 1/2 and 3/2 couplings being more than double the 
PDG estimate. 



TABLE III: Resonance parameters for N* and A* from the SAID fit to the ttN data [51.], helicity 
amplitudes ^1/2 and from the FA08 solution, MAID07 determination [46], and average values 



from Ref. 



Resonance 


ttN said 


A/2 (GeV-1/2 X 10-3) 


A3/2 (GeV-1/2 X 10-3) 




Wr (MeV) 


r (MeV) 


r./r 


FA08 


MAID07 


PDG 


FA08 


MAID07 


PDG 


7V(1535)5ii 


1547 


188 


0.36 


100.9±3.0 


66 


90±30 








A^(1650)5ii 


1635 


115 


1.00 


9.0±9.1 


33 


53±16 








A^(1440)Pii 


1485 


284 


0.79 


-56.4±1.7 


-61 


-65±4 








A^(1720)Pi3 


1764 


210 


0.09 


90.5±3.3 


73 


18±30 


-36.0±3.9 


-11 


-19±20 


iV(1520)Z)i3 


1515 


104 


0.63 


-26±1.5 


-27 


-24±9 


141.2±1.7 


161 


166±5 


iV(1675)Z)i5 


1674 


147 


0.39 


14.9±2.1 


15 


19±8 


18.4±2.1 


22 


15±9 


Af(1680)Fi5 


1680 


128 


0.70 


-17.6±1.5 


-25 


-15±6 


134.2±1.6 


134 


133±12 


A( 1620) 531 


1615 


147 


0.32 


47.2±2.3 


66 


27±11 








A(1232)P33 


1233 


119 


1.00 


-139.6±1.8 


-140 


-135±6 


-258.9±2.3 


-265 


-250±8 


A(1700)Z)33 


1695 


376 


0.16 


118.3±3.3 


226 


104±15 


110.0±3.5 


210 


85±22 


A(1905)F35 


1858 


321 


0.12 


11.4±8.0 


18 


26±11 


-51.0±8.0 


-28 


-45±20 


A(1950)F37 


1921 


271 


0.47 


-71.5±1.8 


-94 


-76±12 


-94.7±1.8 


-121 


-97±10 
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VIII. CONCLUSION 



Differential cross sections for vr"*" meson photopro duct ion on the proton via the reaction 
'J p —>■ n 7T^ have been determined with a tagged-photon beam for incident photon energies 
from 0.725 to 2.875 GeV. All derived cross sections were based on a vr+n missing mass 
reconstruction. The relative cross sections were determined from yields derived from a peak 
isolated above a well-determined background, using Monte Carlo simulations to determine 
the TT"*" acceptance in the CLAS spectrometer. The relative differential cross sections were 
converted to absolute differential cross sections by measurements of the incident photon flux. 

These data have been included in a SAID multipole analysis, resulting in a new SAID 
solution, FA08. Comparisons to earlier SAID fits and a fit from the Mainz group show that 
the new solution is much more satisfactory at higher energies. Although resonance couplings 
have not changed significantly with the addition of these cross sections to the world data 
set, significant changes have occurred in the high-energy behavior of the SAID cross-section 
predictions and amplitudes, as can be seen in Fig. [5] for the cross-section and Fig. [H for 
the single-pion contribution to the GDH sum rule. Further improvement will be possible 
with future measurements of spin observables for the photoproduction process that can be 
expected from FROST |59j and the g8b CLAS running period |60 |. 
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